Substantial viral diversity in bats and rodents from East Africa: insights into evolution, recombination, and cocirculation.

BACKGROUND: Zoonotic viruses cause substantial public health and socioeconomic problems worldwide. Understanding how viruses evolve and spread within and among wildlife species is a critical step when aiming for proactive identification of viral threats to prevent future pandemics. Despite the many proposed factors influencing viral diversity, the genomic diversity and structure of viral communities in East Africa are largely unknown. RESULTS: Using 38.3 Tb of metatranscriptomic data obtained via ultradeep sequencing, we screened vertebrate-associated viromes from 844 bats and 250 rodents from Kenya and Uganda collected from the wild. The 251 vertebrate-associated viral genomes of bats (212) and rodents (39) revealed the vast diversity, host-related variability, and high geographic specificity of viruses in East Africa. Among the surveyed viral families, Coronaviridae and Circoviridae showed low host specificity, high conservation of replication-associated proteins, high divergence among viral entry proteins, and frequent recombination. Despite major dispersal limitations, recurrent mutations, cocirculation, and occasional gene flow contribute to the high local diversity of viral genomes. CONCLUSIONS: The present study not only shows the landscape of bat and rodent viromes in this zoonotic hotspot but also reveals genomic signatures driven by the evolution and dispersal of the viral community, laying solid groundwork for future proactive surveillance of emerging zoonotic pathogens in wildlife. Video Abstract.

Early detection of emerging viral variants through analysis of community structure of coordinated substitution networks.

The emergence of viral variants with altered phenotypes is a public health challenge underscoring the need for advanced evolutionary forecasting methods. Given extensive epistatic interactions within viral genomes and known viral evolutionary history, efficient genomic surveillance necessitates early detection of emerging viral haplotypes rather than commonly targeted single mutations. Haplotype inference, however, is a significantly more challenging problem precluding the use of traditional approaches. Here, using SARS-CoV-2 evolutionary dynamics as a case study, we show that emerging haplotypes with altered transmissibility can be linked to dense communities in coordinated substitution networks, which become discernible significantly earlier than the haplotypes become prevalent. From these insights, we develop a computational framework for inference of viral variants and validate it by successful early detection of known SARS-CoV-2 strains. Our methodology offers greater scalability than phylogenetic lineage tracing and can be applied to any rapidly evolving pathogen with adequate genomic surveillance data.

Genomic analyses of Symbiomonas scintillans show no evidence for endosymbiotic bacteria but does reveal the presence of giant viruses.

Symbiomonas scintillans Guillou et Chrétiennot-Dinet, 1999 is a tiny (1.4 µm) heterotrophic microbial eukaryote. The genus was named based on the presence of endosymbiotic bacteria in its endoplasmic reticulum, however, like most such endosymbionts neither the identity nor functional association with its host were known. We generated both amplification-free shotgun metagenomics and whole genome amplification sequencing data from S. scintillans strains RCC257 and RCC24, but were unable to detect any sequences from known lineages of endosymbiotic bacteria. The absence of endobacteria was further verified with FISH analyses. Instead, numerous contigs in assemblies from both RCC24 and RCC257 were closely related to prasinoviruses infecting the green algae Ostreococcus lucimarinus, Bathycoccus prasinos, and Micromonas pusilla (OlV, BpV, and MpV, respectively). Using the BpV genome as a reference, we assembled a near-complete 190 kbp draft genome encoding all hallmark prasinovirus genes, as well as two additional incomplete assemblies of closely related but distinct viruses from RCC257, and three similar draft viral genomes from RCC24, which we collectively call SsVs. A multi-gene tree showed the three SsV genome types branched within highly supported clades with each of BpV2, OlVs, and MpVs, respectively. Interestingly, transmission electron microscopy also revealed a 190 nm virus-like particle similar the morphology and size of the endosymbiont originally reported in S. scintillans. Overall, we conclude that S. scintillans currently does not harbour an endosymbiotic bacterium, but is associated with giant viruses.

Unveiling hidden structural patterns in the SARS-CoV-2 genome: Computational insights and comparative analysis.

SARS-CoV-2, the causative agent of COVID-19, is known to exhibit secondary structures in its 5' and 3' untranslated regions, along with the frameshifting stimulatory element situated between ORF1a and 1b. To identify additional regions containing conserved structures, we utilized a multiple sequence alignment with related coronaviruses as a starting point. We applied a computational pipeline developed for identifying non-coding RNA elements. Our pipeline employed three different RNA structural prediction approaches. We identified forty genomic regions likely to harbor structures, with ten of them showing three-way consensus substructure predictions among our predictive utilities. We conducted intracomparisons of the predictive utilities within the pipeline and intercomparisons with four previously published SARS-CoV-2 structural datasets. While there was limited agreement on the precise structure, different approaches seemed to converge on regions likely to contain structures in the viral genome. By comparing and combining various computational approaches, we can predict regions most likely to form structures, as well as a probable structure or ensemble of structures. These predictions can be used to guide surveillance, prophylactic measures, or therapeutic efforts. Data and scripts employed in this study may be found at https://doi.org/10.5281/zenodo.8298680.

Conserved structures and dynamics in 5'-proximal regions of Betacoronavirus RNA genomes.

Betacoronaviruses are a genus within the Coronaviridae family of RNA viruses. They are capable of infecting vertebrates and causing epidemics as well as global pandemics in humans. Mitigating the threat posed by Betacoronaviruses requires an understanding of their molecular diversity. The development of novel antivirals hinges on understanding the key regulatory elements within the viral RNA genomes, in particular the 5'-proximal region, which is pivotal for viral protein synthesis. Using a combination of cryo-electron microscopy, atomic force microscopy, chemical probing, and computational modeling, we determined the structures of 5'-proximal regions in RNA genomes of Betacoronaviruses from four subgenera: OC43-CoV, SARS-CoV-2, MERS-CoV, and Rousettus bat-CoV. We obtained cryo-electron microscopy maps and determined atomic-resolution models for the stem-loop-5 (SL5) region at the translation start site and found that despite low sequence similarity and variable length of the helical elements it exhibits a remarkable structural conservation. Atomic force microscopy imaging revealed a common domain organization and a dynamic arrangement of structural elements connected with flexible linkers across all four Betacoronavirus subgenera. Together, these results reveal common features of a critical regulatory region shared between different Betacoronavirus RNA genomes, which may allow targeting of these RNAs by broad-spectrum antiviral therapeutics.

Viral genome sequencing to decipher in-hospital SARS-CoV-2 transmission events.

The SARS-CoV-2 pandemic has highlighted the need to better define in-hospital transmissions, a need that extends to all other common infectious diseases encountered in clinical settings. To evaluate how whole viral genome sequencing can contribute to deciphering nosocomial SARS-CoV-2 transmission 926 SARS-CoV-2 viral genomes from 622 staff members and patients were collected between February 2020 and January 2021 at a university hospital in Munich, Germany, and analysed along with the place of work, duration of hospital stay, and ward transfers. Bioinformatically defined transmission clusters inferred from viral genome sequencing were compared to those inferred from interview-based contact tracing. An additional dataset collected at the same time at another university hospital in the same city was used to account for multiple independent introductions. Clustering analysis of 619 viral genomes generated 19 clusters ranging from 3 to 31 individuals. Sequencing-based transmission clusters showed little overlap with those based on contact tracing data. The viral genomes were significantly more closely related to each other than comparable genomes collected simultaneously at other hospitals in the same city (n = 829), suggesting nosocomial transmission. Longitudinal sampling from individual patients suggested possible cross-infection events during the hospital stay in 19.2% of individuals (14 of 73 individuals). Clustering analysis of SARS-CoV-2 whole genome sequences can reveal cryptic transmission events missed by classical, interview-based contact tracing, helping to decipher in-hospital transmissions. These results, in line with other studies, advocate for viral genome sequencing as a pathogen transmission surveillance tool in hospitals.

Oxygen-dependent histone lysine demethylase 4 restricts hepatitis B virus replication.

Mammalian cells have evolved strategies to regulate gene expression when oxygen is limited. Hypoxia-inducible factors (HIF) are the major transcriptional regulators of host gene expression. We previously reported that HIFs bind and activate hepatitis B virus (HBV) DNA transcription under low oxygen conditions; however, the global cellular response to low oxygen is mediated by a family of oxygenases that work in concert with HIFs. Recent studies have identified a role for chromatin modifiers in sensing cellular oxygen and orchestrating transcriptional responses, but their role in the HBV life cycle is as yet undefined. We demonstrated that histone lysine demethylase 4 (KDM4) can restrict HBV, and pharmacological or oxygen-mediated inhibition of the demethylase increases viral RNAs derived from both episomal and integrated copies of the viral genome. Sequencing studies demonstrated that KDM4 is a major regulator of the hepatic transcriptome, which defines hepatocellular permissivity to HBV infection. We propose a model where HBV exploits cellular oxygen sensors to replicate and persist in the liver. Understanding oxygen-dependent pathways that regulate HBV infection will facilitate the development of physiologically relevant cell-based models that support efficient HBV replication.

Bipartite viral RNA genome heterodimerization influences genome packaging and virion thermostability.

Multi-segmented viruses often multimerize their genomic segments to ensure efficient and stoichiometric packaging of the correct genetic cargo. In the bipartite Nodaviridae family, genome heterodimerization is also observed and conserved among different species. However, the nucleotide composition and biological function for this heterodimer remain unclear. Using Flock House virus as a model system, we developed a next-generation sequencing approach ("XL-ClickSeq") to probe heterodimer site sequences. We identified an intermolecular base-pairing site which contributed to heterodimerization in both wild-type and defective virus particles. Mutagenic disruption of this heterodimer site exhibited significant deficiencies in genome packaging and encapsidation specificity to viral genomic RNAs. Furthermore, the disruption of this intermolecular interaction directly impacts the thermostability of the mature virions. These results demonstrate that the intermolecular RNA-RNA interactions within the encapsidated genome of an RNA virus have an important role on virus particle integrity and thus may impact its transmission to a new host.IMPORTANCEFlock House virus is a member of Nodaviridae family of viruses, which provides a well-studied model virus for non-enveloped RNA virus assembly, cell entry, and replication. The Flock House virus genome consists of two separate RNA molecules, which can form a heterodimer upon heating of virus particles. Although similar RNA dimerization is utilized by other viruses (such as retroviruses) as a packaging mechanism and is conserved among Nodaviruses, the role of heterodimerization in the Nodavirus replication cycle is unclear. In this research, we identified the RNA sequences contributing to Flock House virus genome heterodimerization and discovered that such RNA-RNA interaction plays an essential role in virus packaging efficiency and particle integrity. This provides significant insight into how the interaction of packaged viral RNA may have a broader impact on the structural and functional properties of virus particles.

Functional properties of measles virus proteins derived from a subacute sclerosing panencephalitis patient who received repeated remdesivir treatments.

Subacute sclerosing panencephalitis (SSPE) is a rare but fatal late neurological complication of measles, caused by persistent measles virus (MeV) infection of the central nervous system. There are no drugs approved for the treatment of SSPE. Here, we followed the clinical progression of a 5-year-old SSPE patient after treatment with the nucleoside analog remdesivir, conducted a post-mortem evaluation of the patient's brain, and characterized the MeV detected in the brain. The quality of life of the patient transiently improved after the first two courses of remdesivir, but a third course had no further clinical effect, and the patient eventually succumbed to his condition. Post-mortem evaluation of the brain displayed histopathological changes including loss of neurons and demyelination paired with abundant presence of MeV RNA-positive cells throughout the brain. Next-generation sequencing of RNA isolated from the brain revealed a complete MeV genome with mutations that are typically detected in SSPE, characterized by a hypermutated M gene. Additional mutations were detected in the polymerase (L) gene, which were not associated with resistance to remdesivir. Functional characterization showed that mutations in the F gene led to a hyperfusogenic phenotype predominantly mediated by N465I. Additionally, recombinant wild-type-based MeV with the SSPE-F gene or the F gene with the N465I mutation was no longer lymphotropic but instead efficiently disseminated in neural cultures. Altogether, this case encourages further investigation of remdesivir as a potential treatment of SSPE and highlights the necessity to functionally understand SSPE-causing MeV.IMPORTANCEMeasles virus (MeV) causes acute, systemic disease and remains an important cause of morbidity and mortality in humans. Despite the lack of known entry receptors in the brain, MeV can persistently infect the brain causing the rare but fatal neurological disorder subacute sclerosing panencephalitis (SSPE). SSPE-causing MeVs are characterized by a hypermutated genome and a hyperfusogenic F protein that facilitates the rapid spread of MeV throughout the brain. No treatment against SSPE is available, but the nucleoside analog remdesivir was recently demonstrated to be effective against MeV in vitro. We show that treatment of an SSPE patient with remdesivir led to transient clinical improvement and did not induce viral escape mutants, encouraging the future use of remdesivir in SSPE patients. Functional characterization of the viral proteins sheds light on the shared properties of SSPE-causing MeVs and further contributes to understanding how those viruses cause disease.

Revisiting host identification of bacteriophage Enc34: from biochemical to molecular.

In our 2012 genome announcement (J Virol 86:11403-11404, 2012, https://doi.org/10.1128/JVI.01954-12), we initially identified the host bacterium of bacteriophage Enc34 as Enterobacter cancerogenus using biochemical tests. However, later in-house DNA sequencing revealed that the true host is a strain of Hafnia alvei. Capitalizing on our new DNA-sequencing capabilities, we also refined the genomic termini of Enc34, confirming a 60,496-bp genome with 12-nucleotide 5' cohesive ends. IMPORTANCE: Our correction reflects the evolving landscape of bacterial identification, where molecular methods have supplanted traditional biochemical tests. This case underscores the significance of revisiting past identifications, as seemingly known bacterial strains may yield unexpected discoveries, necessitating essential updates to the scientific record. Despite the host identity correction, our genome announcement retains importance as the first complete genome sequence of a Hafnia alvei bacteriophage.

Cryptic proteins translated from deletion-containing viral genomes dramatically expand the influenza virus proteome.

Productive infections by RNA viruses require faithful replication of the entire genome. Yet many RNA viruses also produce deletion-containing viral genomes (DelVGs), aberrant replication products with large internal deletions. DelVGs interfere with the replication of wild-type virus and their presence in patients is associated with better clinical outcomes. The DelVG RNA itself is hypothesized to confer this interfering activity. DelVGs antagonize replication by out-competing the full-length genome and triggering innate immune responses. Here, we identify an additionally inhibitory mechanism mediated by a new class of viral proteins encoded by DelVGs. We identified hundreds of cryptic viral proteins translated from DelVGs. These DelVG-encoded proteins (DPRs) include canonical viral proteins with large internal deletions, as well as proteins with novel C-termini translated from alternative reading frames. Many DPRs retain functional domains shared with their full-length counterparts, suggesting they may have activity during infection. Mechanistic studies of DPRs derived from the influenza virus protein PB2 showed that they poison replication of wild-type virus by acting as dominant-negative inhibitors of the viral polymerase. These findings reveal that DelVGs have a dual inhibitory mechanism, acting at both the RNA and protein level. They further show that DPRs have the potential to dramatically expand the functional proteomes of diverse RNA viruses.

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing.

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

Dinucleotide biases in the genomes of prokaryotic and eukaryotic dsDNA viruses and their hosts.

The genomes of cellular organisms display CpG and TpA dinucleotide composition biases. Such biases have been poorly investigated in dsDNA viruses. Here, we show that in dsDNA virus, bacterial, and eukaryotic genomes, the representation of TpA and CpG dinucleotides is strongly dependent on genomic G + C content. Thus, the classical observed/expected ratios do not fully capture dinucleotide biases across genomes. Because a larger portion of the variance in TpA frequency was explained by G + C content, we explored which additional factors drive the distribution of CpG dinucleotides. Using the residuals of the linear regressions as a measure of dinucleotide abundance and ancestral state reconstruction across eukaryotic and prokaryotic virus trees, we identified an important role for phylogeny in driving CpG representation. Nonetheless, phylogenetic ANOVA analyses showed that few host associations also account for significant variations. Among eukaryotic viruses, most significant differences were observed between arthropod-infecting viruses and viruses that infect vertebrates or unicellular organisms. However, an effect of viral DNA methylation status (either driven by the host or by viral-encoded methyltransferases) is also likely. Among prokaryotic viruses, cyanobacteria-infecting phages resulted to be significantly CpG-depleted, whereas phages that infect bacteria in the genera Burkolderia and Staphylococcus were CpG-rich. Comparison with bacterial genomes indicated that this effect is largely driven by the general tendency for phages to resemble the host's genomic CpG content. Notably, such tendency is stronger for temperate than for lytic phages. Our data shed light into the processes that shape virus genome composition and inform manipulation strategies for biotechnological applications.

Genetic diversity of RNA viruses infecting invertebrate pests of rice.

Invertebrate species are a natural reservoir of viral genetic diversity, and invertebrate pests are widely distributed in crop fields. However, information on viruses infecting invertebrate pests of crops is limited. In this report, we describe the deep metatranscriptomic sequencing of 88 invertebrate samples covering all major invertebrate pests in rice fields. We identified 296 new RNA viruses and 13 known RNA viruses. These viruses clustered within 31 families, with many highly divergent viruses constituting potentially new families and genera. Of the identified viruses, 13 RNA viruses clustered within the Fiersviridae family of bacteriophages, and 48 RNA viruses clustered within families and genera of mycoviruses. We detected known rice viruses in novel invertebrate hosts at high abundances. Furthermore, some novel RNA viruses have genome structures closely matching to known plant viruses and clustered within genera of several plant virus species. Forty-five potential insect pathogenic RNA viruses were detected in invertebrate species. Our analysis revealed that host taxonomy plays a major role and geographical location plays an important role in structuring viral diversity. Cross-species transmission of RNA viruses was detected between invertebrate hosts. Newly identified viral genomes showed extensive variation for invertebrate viral families or genera. Together, the large-scale metatranscriptomic analysis greatly expands our understanding of RNA viruses in rice invertebrate species, the results provide valuable information for developing efficient strategies to manage insect pests and virus-mediated crop diseases.

PhageScope: a well-annotated bacteriophage database with automatic analyses and visualizations.

Bacteriophages are viruses that infect bacteria or archaea. Understanding the diverse and intricate genomic architectures of phages is essential to study microbial ecosystems and develop phage therapy strategies. However, the existing phage databases are short of meticulous annotations. To this end, we propose PhageScope (https://phagescope.deepomics.org), an online phage database with comprehensive annotations. PhageScope harbors a collection of 873 718 phage sequences from various sources. Applying fifteen state-of-the-art tools to perform systematic annotations and analyses, PhageScope provides annotations on genome completeness, host range, lifestyle information, taxonomy classification, nine types of structural and functional genetic elements, and three types of comparative genomic studies for curated phages. Additionally, PhageScope incorporates automatic analyses and visualizations for curated and customized phages, serving as an efficient platform for phage study.

Structural similarity of human papillomavirus E4 and polyomaviral VP4 exhibited by genomic analysis of the common kestrel (Falco tinnunculus) polyomavirus.

Polyomaviruses are widely distributed viruses of birds that may induce developmental deformities and internal organ disorders primarily in nestlings. In this study, polyomavirus sequence was detected in kidney and liver samples of a common kestrel (Falco tinnunculus) that succumbed at a rescue station in Hungary. The amplified 5025 nucleotide (nt) long genome contained the early (large and small T antigen, LTA and STA) and late (viral proteins, VP1, VP2, VP3) open reading frames (ORFs) typical for polyomaviruses. One of the additional putative ORFs (named VP4) showed identical localization with the VP4 and ORF-X of gammapolyomaviruses, but putative splicing sites could not be found in its sequence. Interestingly, the predicted 123 amino acid (aa) long protein sequence showed the highest similarity with human papillomavirus E4 early proteins in respect of the aa distribution and motif arrangement implying similar functions. The LTA of the kestrel polyomavirus shared <59.2% nt and aa pairwise identity with the LTA sequence of other polyomaviruses and formed a separated branch in the phylogenetic tree among gammapolyomaviruses. Accordingly, the kestrel polyomavirus may be the first member of a novel species within the Gammapolyomavirus genus, tentatively named Gammapolyomavirus faltin.

The most exposed regions of SARS-CoV-2 structural proteins are subject to strong positive selection and gene overlap may locally modify this behavior.

The SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) pandemic that emerged in 2019 has been an unprecedented event in international science, as it has been possible to sequence millions of genomes, tracking their evolution very closely. This has enabled various types of secondary analyses of these genomes, including the measurement of their sequence selection pressure. In this work, we have been able to measure the selective pressure of all the described SARS-CoV-2 genes, even analyzed by sequence regions, and we show how this type of analysis allows us to separate the genes between those subject to positive selection (usually those that code for surface proteins or those exposed to the host immune system) and those subject to negative selection because they require greater conservation of their structure and function. We have also seen that when another gene with an overlapping reading frame appears within a gene sequence, the overlapping sequence between the two genes evolves under a stronger purifying selection than the average of the non-overlapping regions of the main gene. We propose this type of analysis as a useful tool for locating and analyzing all the genes of a viral genome when an adequate number of sequences are available.IMPORTANCEWe have analyzed the selection pressure of all severe acute respiratory syndrome coronavirus 2 genes by means of the nonsynonymous (Ka) to synonymous (Ks) substitution rate. We found that protein-coding genes are exposed to strong positive selection, especially in the regions of interaction with other molecules (host receptor and genome of the virus itself). However, overlapping coding regions are more protected and show negative selection. This suggests that this measure could be used to study viral gene function as well as overlapping genes.

A virus from Aspergillus cibarius with features of alpha- and betachrysoviruses.

From the ascomycete Aspergillus cibarius strain NW-FVA 2590, which was originally isolated from a root, associated with stem collar necrosis of Fraxinus excelsior L., a novel virus was isolated and characterized. Its genome is encoded on three monocistronic dsRNA segments ranging from 3683 bp (dsRNA 1) over 3093 (dsRNA 2) to 2902 bp (dsRNA 3), which are packed in isometric particles of around 35 nm. While the viral RdRp (P1) is encoded on segment 1, protein sequencing showed that two more structural proteins are present which are translated from dsRNA 2 (P2) and dsRNA 3 (P3) and possibly form the viral capsid. Additionally, P2 and P3 may undergo posttranslational modifications since the detected proteins bands deviated from the calculated sizes. Due to its phylogenetic position, the novel virus was grouped in the family of Chrysoviridae and was tentatively denominated as Aspergillus cibarius chrysovirus 1 (AcCV1). Due to its composition, biological properties and phylogenetic position, distant from the genera Alphachrysovirus and Betachrysovirus, we suggest to position AcCV1 in a proposed genus "Gammachrysovirus".

First report of the whole-genome sequence analysis of avian rotavirus A from Japanese chickens.

Rotavirus A infects many mammalian species, including humans and causes diarrhea and gastrointestinal diseases. The virus also infects various bird species, including chickens, although information of avian rotavirus A (ARVA) infection in chicken populations in Japan is scarce. In this study, we report for the first time the whole-genome sequences of ARVA strains from Japanese chicken populations. The virus strains were inoculated to MA104 cells and cultured viruses were used to obtain the sequences with the MiSeq system, and genetic analysis demonstrated the genotype constellation of G19-P[30]-I11-R6-C6-M7-A16-N6-T8-E10-H8 of the Japanese chicken ARVA isolates. Phylogenetic analyses demonstrated that the VP1, VP2, VP3, VP4, VP7, NSP2, and NSP4 coding gene sequences of the Japanese strains were closer to those of Korean than the European ARVA strains, although such relationship was not clear for other genes. The data suggest that the Japanese ARVA strains and the ones in Korea have genetically close relationship, although the origin is not clear at this point. Further information including the whole-genome sequences of the Korean strains and sequences of other Japanese chicken ARVA strains will be necessary for elucidation of their origin.